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Abstract. Using adaptive optics at the ESO 3.6m telescope, we 
obtained diffraction limited JHK-images of the region around 
the Herbig AeBe star NX Pup. We clearly resolved the close 
companion (sep. 0'.'128) to NX Pup - originally discovered by 
HST - and measured its JHK magnitudes. A third object at a 
separation of 7'.'0 from NX Pup was identified as a classical T 
Tauri star so that NX Pup may in fact form a hierarchical triple 
system. We discuss the evolutionary status of these stars and 
derive estimates for their spectral types, luminosities, masses 
and ages. 
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1. Introduction 

Cometary Globules (CGs) are cool clouds with compact, dusty, 
opaque heads, and faint tails stretching away from the head 
(e.g., Harju et al. 1990). They were first described and defined 
as a group by Hawarden and Brand (1976) based on Schmidt 
plates of the Gum Nebula and NGC 5367. Many CGs are found 
in the Gum Nebula, a large region of ionized gas at a distance of 
approximately 450pc. In this nebula the isolated neutral glob- 
ules clearly contrast with the surrounding hot ionized media 
(Brand et al. 1983). 



Send offprint requests to: W. Brandner 

* based on observations obtained at the European Southern Obser- 
vatory, La Silla, Chile 



The heads of most CGs point to an apparent center in the 
Gum Nebula, suggesting a common origin of the cometary 
shape of these globules (Zealey et al. 1982). Different mech- 
anisms have been proposed to explain the origin of the tail like 
structure of the CGs, such as interactions of dense cloud cores 
with shock waves from a SN explosion (Brand et al. 1983) 
or ionization shock fronts from nearby OB stars (Reipurth 
1983). IRAS observations show an enhanced star formation 
rate in the CGs of the Gum Nebula with respect to compar- 
ison fields (Bhatt 1993). Star formation may have been in- 
duced through shocks triggered by the same event that ion- 
ized the Gum Nebula complex. Reipurth and Pettersson (1993) 
identified eight late-type Ha emission stars in association with 
CG4/CG6/Sal01 and CGI 3. They argue that star formation has 
occurred in the Gum Nebula for more than 10 6 years. A map- 
ping of CGI in molecular lines has been carried out by Harju 
et al. (1990). They find that CGI is not in dynamical equilib- 
rium and in CO shows a high (5km/s) velocity component near 
its head. They conclude that this is shocked material related to 
star forming processes. 

Located at the edge of CGI is the Herbig AeBe star NX Pup 
(Irvine 1975). NX Pup is a variable star (Hoffmeister 1949, 
Strohmeier et al. 1964). It shows Ha in emission (Wackerling 
1970) and has a strong UV excess (de Boer 1977) and IR ex- 
cess (Brand et al. 1983, Reipurth 1983). Its evolutionary sta- 
tus is controversial. Brand et al. and Reipurth find an MK type 
of F0-2III, a total luminosity of 40-50 L©, and an age of 10 6 
years, whereas (on the basis of UV observations from space) de 
Boer (1977) and Tjin A Djie et al. (1984) assign an MK type of 
AO-1 III. Tjin A Djie et al. and The and Molster (1994) derive 
a total luminosity between 100 and 160L Q and an age of 8 TO 5 
years for NX Pup. The age of NX Pup agrees very well with the 
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age of CGI derived from dynamical studies (e.g., Harju et al. 
1990). 

Observations with the Fine Guidance Sensor system aboard 
the Hubble Space Telescope reveal that NX Pup (HIC 35488) 
is a close visual binary with a separation of 0'.'126 (Bernacca et 
al. 1993). Thus all previous evolutionary interpretations should 
be reevaluated as the total observed luminosity in fact comes 
from two stars. 

Using the ESO adaptive optics (AO) imaging system 
COME-ON+ (CO+) in combination with the SHARP camera 
we started a diffraction limited imaging survey of intermediate- 
mass pre-main sequence (PMS) stars listed in "A new catalogue 
of members and candidate members of the Herbig Ae/Be stel- 
lar group" (The et al. 1994) with the aim to search for close IR 
companions. During the first part of this project we obtained 
JHK images for the NX Pup system. The CLEANed JHK im- 
ages clearly resolve the close companion, here referred to as 
NX Pup B, detected in the optical by the HST 

We also found a third, faint PMS star at a distance of 7'.'0 
from NX Pup A/B, which will be referred to as NX Pup C in 
the following. 



Table 1. Journal of observations 



telescope/instr. 


date 


filter/wavelength 


ESO lm/vis. phot. 


1.-3.12.1992 


V 


SAT (LTPV)/vis. phot. 


12.92-2.93 


y 


ES0 3.6m/CO+ 


1.1.1994 


J,H,K 


Danish 1.5m/CCD-camera 


8.1.1994 


U,B,V,R,i,Ha 


ESO 1.5m/B&C 


26.1.1994 


380-760nm 


NTT/EMMI 


20.3.1994 


600-725nm 




1 I Lo.o 

-0.4 -0.2 -0.0 0.2 0.4 
arcsec 



Fig. 1. Deconvolved JHK images of NX Pup obtained with the adap- 
tive optics system CO+ at the ESO 3.6m telescope in 1994 Jan 1. The 
images were rebinned by a factor of 4 (see text for more details). The 
components A and B (sep. 07128, PA 62.4°) are clearly resolved. At 
a distance of 450pc 0'.'128 correspond to a projected separation of 
58AU. The faint feature in the J image north-west of NX Pup A is 
an artefact of the image deconvolution. A logarithmic gray scale was 
used. Component C is outside these frames. North is up and east is to 
the left. 



technique, i.e., coadding the individual images after shifting 
them so that the centre of the point spread functions coincide. 



2. Observations and results 

2.1. Near-infrared Imaging 

JHK images of NX Pup with high spatial resolution were ob- 
tained in 1994 Jan 1 at the ESO 3.6m telescope. We used 
the COME-ON+ adaptive optics system (Beuzit and Rousset, 
1994) in combination with the SHARP II (System for High 
Angular Resolution Pictures, see Hofmann et al. 1992) cam- 
era from the Max Planck Institute for Extraterrestrial Physics 
(MPE). The SHARP camera is equipped with a Rockwell 
NICMOS-3 array. The image scale was O'.'050/pixel, giving a 
field size of 13" x 13". The total exposure time on NX Pup 
was 300s in each JHK filter, and similar exposures were ob- 
tained on the nearby sky for sky subtraction. Immediately after 
NX Pup, we observed a reference point source 10' away (star 
no. 985 in the HST Guide Star Catalogue) which was later used 
for image deconvolution. The photometric standard HR 3421 
was observed under the same conditions to get absolute JHK 
photometry. 

The individual exposures in each filter, 5 in J, 150 in H, 
and 600 in K, were coadded using the "shift and add" (SAA) 



Table 2. Gain in spatial resolution (FWHM and Strehl ratio, SR) 
through shift and add (SAA) and SAA with image selection (IS). 





J 


H 


K 


FWHM 


07214 


0'.'172 


0'.'156 


FWHM (SAA) 


07209 


07157 


07147 


FWHM (SAA+IS) 


a 


07138 


07143 


SR 


4% 


8% 


25% 


SR (SAA) 


4% 


9% 


27% 


SR (SAA+IS) 




12% 


31% 


selection rate (IS) 




20% 


40% 



a no IS (only 5 individual exposures) 



All AO systems do only partially compensate the wavefront 
distorsions. "The performance of an AO system is best defined 
in terms of the maximum intensity obtained in a point source 
image. The ratio of this intensity to the maximum intensity in 
the diffraction-limited image (Airy pattern) is called the Strehl 
ratio" (Roddier 1992). The uncompensated power forms a halo 
around each object. 
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If a sufficient number (> 50) of individual exposures is 
available, one can gain in spatial resolution by selecting images 
based on a Strehl ratio (SR) criteria. Furthermore, if the single 
frame exposure times are short, the residual tilt between the 
individual images is reduced by applying an SAA algorithm. 
While the latter already gives some improvements mostly in 
the FWHM, the results of image selection are more significant 
with respect to the Strehl ratio. However, the rejection rate for 
image selection must be limited in order to keep a sufficient 
signal to noise ratio. Table 2 shows the improvement in reso- 
lution through SAA and image selection. Out of 600 K band 
frames 240 were selected and coadded, improving the Strehl 
ratio from 25% to 31%. Selecting 30 images out of the 150 H 
band exposures improved the Strehl ratio to 12% - compared 
to 8% without image selection. 

All these operations were performed with the "local" IRAF 
package cl28 developed by E. Tessier at the Observatory of 
Grenoble 1 . See also Tessier et al. (1995) for a more detailed 
description of the data reduction . 

On the shift-and-add image in the K band, NX PUP C itself 
appears to be elongated in the east-west direction (FWHM in 
east-west direction and north-south direction are 0'.'27 x 0'.'18, 
compared to 0'.'16 x 0'.'16 for the reference source). This elon- 
gation is less pronounced in the H band, and is not seen in the 
J band. 

NX Pup C is 7" away from NX Pup AB, which was used 
for wave-front sensing. Therefore, we could expect effects from 
angular anisoplanatism (Wilson and Jenkins 1995). Because of 
this and the lower signal to noise ratio of NX Pup C, we did not 
try to deconvolve it. 

In the following we will discuss briefly two kinds of aniso- 
planatism - angular and temporal - which cause image elon- 
gations. One effect of angular anisoplanatism is that the Strehl 
ratio drops and the PSF gets wider with increasing angular dis- 
tance to the guide star. The FWHM of NX Pup C is systemat- 
ically wider than the FWHM of the calibration source in each 
JHK filter. This might be because of angular anisoplanatism. 
Furthermore, angular anisoplanatism would lead to an elonga- 
tion of the off-axis PSF in the direction of the guide source 
(McClure et al. 1991). In the case of NX Pup AB as reference 
source and NX Pup C as programme star the direction of the 
elongation would be ~ 45°, i.e., from north-east to south-west 
but not in east-west direction. 

Elongations in another direction can be explained by tem- 
poral anisoplanatism due to the wind in the dominant turbulent 
layer (see Roddier et al. 1993 and Wilson and Jenkins 1994). 
The prevalent wind direction on La Silla near the ground is 
north-south, but we do not know the wind direction in the domi- 
nant turbulent layer during our observations. While an off-axis 
PSF should suffer from both effects, the on-axis PSF should 
show pure temporal anisoplanatism. However, the on-axis PSF 
of the reference source does not show any significant devia- 
tion from circular symmetry. Summarizing, neither angular nor 

1 The IRAF package cl28 is available via anonymous ftp from 
hplyot.obspm.fr. in the directory /iraf_hra. 



temporal anisoplanatism are very likely sources of elongation 
and the elongation may therefore come from the object NX 
Pup C itself. To know for sure NX Pup C would have to be 
re-observed. 

The resulting coadded images were then deconvolved using 
a CLEAN algorithm and the nearby reference star as a PSF. The 
central 0'.'5 x 0'.'5 part of the CLEANed JHK images is shown in 
Fig. 1. The close visual companion detected by HST is clearly 
seen in all filters. This is a clear demonstration of the diffraction 
limited imaging capabilities of the CO+ adaptive optics system. 

Astrometry has been performed on the CLEANed images, 
whereas photometry of the various components of the system 
has been performed on the shift-and-add images. We did not 
carry out the photometry on the CLEANed images as decon- 
volution algorithms usually do not conserve the flux distri- 
bution. Absolute JHK photometry was obtained for NX Pup 
A+B and for NX Pup C using the IRAF/APPHOT package, 
while differential photometry between NX Pup A and B was 
performed within the IRAF/DAOPHOT package using IRAF 
scripts kindly provided by A. Tokovinine. The JHK magnitudes 
derived in this manner for NX Pup A, B, and C are listed in Ta- 
ble 3. The overall photometric errors are about 10%. 



Table 3. JHK photometry (1.1.1994, ESO 3.6m/CO+). 
For comparison we also show the V magnitudes. 



filter 


NX Pup A 


NX Pup B 


NX Pup C 


J 


8 m 58 


9^56 


ll m 71 


H 


7 m 43 


8 m 37 


10" 66 


K 


6 m 15 

1 mil oa 


7 m 90 

i /-aiirva 


ictio 



a assuming AV=(T64 (Bernacca et al. 1993; 
1.1.1993) and V=9 m 8 for NX Pup A+B 
6 this paper (8.1.1994) 



2.2. Visual Imaging 

CCD images in the Bessel UBVR, Gunn i, and Ha filters of 
the region around NX Pup were obtained in 1994 Jan 8 with 
the CCD-camera at the Danish 1 .5m telescope at La Silla. This 
camera is equipped with a IK TEK CCD (ESO #28) and has a 
scale of 0'.'38/pixel. 

We obtained two sets of exposures in U,B,V,R,i,Ha - short 
exposures (10s, 5s, Is, Is, Is, 60s) in order not to saturate 
NX Pup AB and long exposures (180s, 20s, 10s, 10s, 10s, 600s) 
in order to get reliable photometry of NX Pup C. The seeing 
was between l'.'O and l'/5. The finding chart for NX Pup AB 
and C is shown in Fig. 2. 

Two standard fields (SA98 and RU149) of the Landolt list 
(Landolt 1992) were observed in order to calibrate the photo- 
metric data. 

The CCD camera at the Danish 1.5m telescope has a shut- 
ter delay time of 30ms. For short exposures of lsec this re- 
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Fig. 2. The PMS stars NX Pup A+B (center) and C in a 10s R exposure 
with the CCD camera at the Danish 1.5m telescope on La Silla. The 
components AB (unresolved) and C are marked. All other stars within 
the field show no excess in Ha and therefore are likely to be field stars. 
North is up and east is to the left. 

suits in a 3% gradient in intensity from the center to the edge 
of the field. All short exposures were corrected for the shutter 
delay. Instrumental magnitudes and transformation equations 
were computed and solved within the IRAF/DIGIPHOT pack- 
age. The overall photometric errors are between 2% and 5%. 
Table 4 lists the results of the visual photometry. 

As NX Pup AB is variable we have to be careful when com- 
bining the observations of different epochs. We do not have 
an absolute V magnitude for the HST observation of NX Pup, 
but one of us (DdW) observed NX Pup with the visual pho- 
tometer at the ESO lm telescope between 1992 Dec 1 and 3. 
Furthermore, NX Pup is one of the target stars in the LTPV 



Table 4. Visual Photometry (8.1.1994, Danish 1.54m) 



filter 


NX Pup AB 


NX Pup C 


U 


10™ 40 


17^52 


B 


KT30 


17 m 46 


V 


9™78 




R 


9 m 46 


14 m 96 


I 


9 m 01 


13 m 68 


R-Ha a 


(T46 


(T68 



a see definition of Ha magnitude in the text 



programme and was observed between 1992 Dec 4 and 22 and 
1993 Jan 21 and Feb 8. PS. The was so kind to provide us 
Str0mgren y band measurements of NX Pup - obtained with 
the Danish 50cm Str0mgren Automatic Telescope (SAT) on La 
Silla, Chile - prior to their publication (Sterken et al. 1995, in 
prep.). NX Pup was getting fainter from December to Febru- 
ary. If we assume the y band measurements to be equivalent 
to the V magnitudes, simple linear interpolation gives us V = 
9™ 8 ± (P2 at the time of the HST observations. Thus NX Pup 
had about the same V magnitude during the CO+ observations 
as during the HST observations. 



2.3. Spectroscopy 

We used the ESO 1.5m and the Boiler and Chivens (B&C) 
spectrograph in 1994 Jan 26 to get long slit spectra of 
NXPupAB and C simultaneously. The B&C spectrograph is 
equipped with a 2K FORD chip (ESO #24) and can be rotated 
as a whole in order to align the slit along a user defined po- 
sition angle. We used grating #23 with an useful wavelength 
range from 380 to 760nm and a sampling of 0. 19nm/pixel. The 
FWHM along the slit of the 2D spectra is 4'.' 3. This relatively 
large value presumably cannot be attributed to the seeing but 
to guiding errors during the 30min exposure. Accordingly, the 
shape of the one dimensional point spread function is better 
described by a Gaussian than by a function resembling a King 
or Moffat profile. Because of the large FWHM and the large 
brightness difference between NXPupAB and C (5.5 mag in 
R), the C component is heavily blended. 

To extract the spectrum of NX Pup C we fit two Gaussians 
simultaneously along each scan in spatial direction of our two 
dimensional spectra. FWHM and separation of the two Gaus- 
sians remain fixed, allowing only their heights and the absolute 
position to vary from scan to scan in the dispersion direction. 
A comparison to the NTT spectrum (see below) shows good 
agreement of the spectrum of NX Pup C in the overlapping re- 
gion, giving us some confidence in our method. 

Spectra with a better spectral resolution (0.06nm / pixel) 
and spatial resolution (FWHM l'/O) were obtained in 1994 Mar 
20 with the ESO Multi-Mode Instrument (EMMI) attached to 
the NTT. We used the red arm of EMMI, which has a 2K TEK 
CCD (ESO #36) and a pixel scale of 0'.'27/pixel. The exposure 
time was 15min and the wavelength range 600 to 725 nm. Be- 
cause of the large brightness difference between NX Pup C and 
AB, at least parts of the spectra of NX Pup AB were saturated 
in both exposures. The saturated regions were excluded from 
the fits and the final extraction. Two spectrophotometric stan- 
dard stars were observed in order to derive the sensitivity curve 
and to achieve a flux calibration. All spectra were extinction 
corrected applying the default La Silla extinction curve. 

Figure 3 shows the medium resolution spectrum of NX Pup 
C near the Ha line. The strong Ha emission and the Li absorp- 
tion can be seen. 
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Fig. 3. Spectrum of NX Pup C, obtained in a 15min exposure with NTT/EMMI in 1994 Mar 20. Several emission (01, Ha) and absorption 
(TiO bands, Ca I) features are marked. Note the strong Li I 670.8nm absorption, which is a sign of youth. 



3. The evolutionary status of NX Pup A+B and C 

3.1. NXPupA+B 

The position angle (PA) between NX Pup AB and C measured 
on the CO+ frames and the CCD frames is 45?3 ± 0?2 and the 
separation 6'.'98±0'.'04. At a distance of 450 pc this corresponds 
to a projected separation of 3140AU, so that NX Pup AB and 
C are unlikely to form a gravitational bound system. On the 
other hand, this resembles what we would expect to see for a 
hierarchical triple system: two stars in a close orbit around each 
other and a third, less massive star in a wide orbit. Only radial 
velocity measurements and proper motion studies could tell us 
whether or not the three stars form a physical triple system. 

NX Pup A and B themselves are separated by 0'.'128 ± 
0'.'008, as estimated from the JHK images, i.e. 58AU pro- 
jected separation. They are almost certainly bound, as has been 
pointed out by Bernacca et al. (1993). The separation and the 
PA of 62?4 ± 5? 7 are in good agreement with the values deter- 
mined by Bernacca et al. (07126 ±07007, 63?4± 1?0). As there 
is only a time span of one year between the HST measurements 
and our CO+ observations, we would not expect signs of orbital 
motion. 

To transform the Ha magnitudes, our Bessell R filter can 
be used as continuum filter. The ratio of our instrumental R and 
Ha filter passbands corresponds to a difference of 3™9 (Grebel 
et al. 1994) in the sense ofR inst ~Ha inst — 3™9 (for stars with- 



out Ha emission). A comparison of the (Ri„ s t — Ha insf + 3™9) 
colours shows that both NX Pup AB and C have a strong ex- 
cess emission in Ha (cf. Table 4). No other stars with R down 
to 19™3 within a 2.'5 x 2f 5 field centered on NX Pup possess a 
significant excess in Ha. 

In H-K vs V-K colour-colour diagrams, NX Pup A and B 
as well as NX Pup C lie outside of the region where reddened 
main sequence stars with normal extinction would be located. 
In the J-H vs H-K diagram NX Pup A falls into the region 
where Herbig AeBe stars can be found. NX Pup B and C lie 
in the region where reddened main sequence stars and T Tauri 
stars with low accretion rates and/or face-on disks (small IR 
excess) are located (Lada and Adams 1992, Strom et al. 1993). 
This indicates that both NX Pup A and B have intrinsic IR ex- 
cesses which may arise from a circumstellar disk or envelope. 
Because of the small separation between A and B, circumstellar 
material either in form of envelopes or disks should be strongly 
disturbed. This has already been pointed out by Henning et 
al. (1994) and could explain the non-detection of NX Pup in 
their 1.3mm continuum search. The spectral energy distribu- 
tion (SED) is plotted in Fig. 4. The major part of the system's 
IR excess emission - already noted by Brand et al. (1983) and 
Reipurth (1983) and fitted by them with a blackbody of T e // 
around 1300K - comes from NX Pup A, as evident in our Fig. 
4. 
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As the reddening vector is almost parallel to the sequence 
of dwarfs and giants with normal colours in the H-K vs V-K 
diagram, we can constrain the spectral type of NX Pup B. If 
extinction along the line of sight is negligible, we derive a spec- 
tral type earlier than K5 for NX Pup B and earlier than G6 if 
A y =l m 5 as for NX Pup A (Hillenbrand et al. 1992). The spec- 
trum obtained with the NTT shows the Ha line of NX Pup AB 
in emission with a central absorption dip. Other lines in emis- 
sion are the Hel line at 587.5nm, the Nal lines at 589.0 and 
589. 6nm, OI at 630. Onm and H/3. Because of the strong satu- 
ration of the CCD in the central part of the spectrum obtained 
at the NTT, we did not try to identify photospheric absorption 
features. The spectrum obtained at the ESO 1.5m shows no ev- 
idence for Li I. Martin (1994) has estimated that, e.g., in the 
combined spectrum of an Herbig AeBe star of spectral type A2 
(50L Q ) and an 1.5M T Tauri type star (10L Q , 10 6 yrs) the Li 
I doublet should be visible with an equivalent width of 0.02nm. 



Table 5. Equivalent width (in nm) of emission (< 0) and absorp- 
tion (> 0) lines in the spectra of NXPupAB and C (28.1.1994, 
ESO 1.5m/B&C; 20.3.1994, NTT/EMMI ) 



-11 



line 


NX Pup AB 


NX Pup C 


Ha 


a 


-2.85 


OI 630.0 


-0.068 


-0.11 


H/3 


0.25 6 


-0.58 


H 7 


0.53 


-0.17 


Li I 670.8 




0.054 



double peaked 
emission core 
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The spectral type of NX Pup A has been under debate. Ob- 
servers in the UV found systematic earlier spectral types (A0- 
A2, de Boer 1977, Tjin A Djie et al. 1984) than observers in 
the visual (F0-F2, Brand et al. 1983, Reipurth 1983). Recently, 
Blondel and Tjin A Djie (1994) suggested that the low resolu- 
tion IUE spectrum of NX Pup A/B can be decomposed into the 
spectrum of a hot boundary layer - situated between accretion 
disk and star - and the photospheric spectrum of an F2 type 
star. 

We use the spectral type - effective temperature calibration 
from de Jager and Nieuwenhuijzen (1987) to determine T e ff. 
If we take the extinction values from Blondel and Tjin A Djie 
(Ay = m 4-0 m 7), assume a spectral type A7-F2 and V=l(ri- 
HT5 for NX Pup A, we get a luminosity of 15-29 L . With 
the knowledge of the luminosity and the effective temperature 
we are able to place NX Pup A on an HR diagram and to de- 
termine its evolutionary status. New sets of pre-main sequence 
evolution tracks have been computed by D' Antona and Mazz- 
itelli (1994). From their tracks based on opacities by Alexander 
et al. (1989) and the convection model from Canuto and Mazz- 
itelli (1990) we get a mass of 2M and an age around 5 x 10 6 
yr (see Table 6 and Fig. 5). 



Fig. 4. Spectral energy distribution AF A of NX Pup A and B (VJHK) 
and C (UB VRIJHK). For comparison we also show the spectral energy 
distribution of a Ml type star of the same apparent V magnitude as NX 
Pup C. Note that the SED is rising towards longer wavelength for NX 
Pup A. The SED of NX Pup B and C peaks near 1.5^im. NX Pup C 
show IR excess and UV excess. The errors in flux are 5% or less. 



The evolutionary status of NX Pup B is less well con- 
strained. From its position in the colour-colour-diagrams and 
the fact, that it is fainter than NX Pup A, a spectral type be- 
tween mid F and late G seems to be likely. Assuming V=10 1 ? 1 7- 
1 1™1 and the same extinction values as for NX Pup A, we get 
a luminosity between 9 and 18 L Q . The dashed box in Fig. 5 
shows that an age from 0.3 to 5 x 10 6 yrs and a mass between 
1.2 and 2.5M Q is possible. 

Spatially resolved observations of NX Pup A and B red- 
ward of 2.2^m would be desirable. However, because of the 
diffraction limit they have to wait for the next generation of 
(southern) telescopes with their 6 to 8 m mirrors. In the mean- 
time speckle observations in the R and I band and speckle spec- 
troscopy could provide us with better estimates on luminosity 
and spectral type of NX Pup A and B. 
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Fig. 5. HR-diagram. The pre-main sequence evolutionary tracks are from D'Antona and Mazzitelli (1994). The solid lines mark the isochrones 
and the zero-age main sequence, the dotted lines are the evolutionary tracks for stars in the mass range from 0. 1 to 2.5Mq . The positions of NX 
Pup A, B, and C in this diagram are marked by boxes. For NX Pup A we assume a spectral type A7-F2, V^O? 1 1-10^5, and A v = (T 1 4-0^7 
(Blondel and Tjin A Djie 1994). For NX Pup B we assume the same extinction, a spectral type F5-G8, and V=l(X n 7-l 1™ 1. The evolutionary 
status of NX Pup C is better defined, yielding an age around 5 x 10 5 yrs and a mass of O.3OM . 



3.2. NX Pup C 

NX Pup C itself is unresolved on the shift-and-add images. In 
the following we regard it as a single star. 

In the spectra of NX Pup C the Li 1 670. 8nm absorption line 
is present with an equivalent width of 0.054nm (Fig. 3). Since 
Li is depleted in the deeper zones of the convection layer of 
low-mass stars, the presence of this line is a sign of youth and 
gives a strong hint that NX Pup C is still in its PMS evolution 
phase. 

The equivalent width of the Ha emission is 2.85nm (Ta- 
ble 5), thus NX Pup C belongs to the group of classical T Tauri 
stars (CTTS). TiO bands at 665.2 and 668.1 nm (typical for M- 
type stars) are visible, whereas TiO bands at 559.8, 562.9 and 
556. 1 nm (present from M2 on) are not visible. The spectral A 



index, which estimates the intensity of the CaH band around 
697.5nm (Kirkpatrick et al. 1991), yields a possible MK type 
between M0.5V and M4III. Taking uncertainties in the extinc- 
tion correction (see below) into account as well as the fact that 
NX Pup C as a PMS star is still in its contraction phase, we get 
a spectral type between M0. 5 and Ml. 5. 

In order to derive NX Pup C's luminosity, an estimate of 
the extinction on the line of sight has to be derived. This is 
not straightforward, since we have three constituents to the ex- 
tinction: foreground extinction, extinction within CGI (radio 
maps from Harju et al. show that the molecular cloud extends 
well beyond NX Pup) and possibly circumstellar extinction. 
Another difficulty when estimating the star's luminosity arises 
from the presence of both an UV excess (see Table 4) and a 
near-IR excess (as revealed by the location of NX Pup C in the 



8 



Brandner et al.: Diffraction limited search for PMS stars near NX Pup 



V-K v. H-K colour-colour diagram, see above). The strong UV 
excess may be explained by the presence of the Balmer jump, 
a feature commonly observed among CTTS (e.g. Valenti et al. 
1993), while the near-IR excess implies that part of the stellar 
radiation is reprocessed within the disk. These excesses and the 
presence of strong Balmer lines indicates ongoing accretion in 
the disk of NX Pup C. Therefore, the observed total luminos- 
ity is the sum of the photospheric luminosity and the accretion 
luminosity. 

As the maximum contribution of the stellar photosphere to 
the total spectral energy distribution (SED) is at about l^m 
(Bertout et al. 1988, Hartigan et al. 1992), the I or J magnitudes 
are the best measures of the stellar photospheric flux. The pres- 
ence of photospheric absorption lines tells us that the veiling 
is not too strong. However, high-resolution spectra would be 
necessary for an exact evaluation of the veiling. 

We use the R-I colour and a spectral type of M1V for 
NX Pup C to compute the extinction. Colours and bolometric 
corrections for M dwarfs are from Hartigan et al. (1994). As in- 
terstellar extinction relations are colour-dependent we use the 
new values for differential and total extinction from Grebel and 
Roberts (1995). Assuming a normal interstellar extinction law 
towards NX Pup and a distance of 450pc, we derive a total ex- 
tinction A/=0.93 (Ay=1.47) and a luminosity of 0.85L Q for 
NX Pup C. If we use the V-R colour instead of R-I, we get 
A/=0.43 (Ay =0.68) and a luminosity of0.45L Q for NX Pup C. 
This illustrates the difficulty in deriving extinction values for 
stars with abnormal spectral energy distributions like classical 
T Tauri stars. 

The spectral energy distribution of NX Pup C is shown in 
Fig. 4 and compared to that of an Ml dwarf. The UV and near- 
IR excess above the photospheric SED are clearly seen, the lat- 
ter amounting to A(K) w 0.5dex, typical for classical T Tauri 
stars (Strom et al. 1989). 

We use again the relation from de Jager and Nieuwenhui- 
jzen (1987) to determine T e // for NX Pup C. In the region be- 
tween MO. 5 and Ml. 5, T e ff is almost independent of the lu- 
minosity class. Taking the uncertainties in the spectral classi- 
fication into account, we get an effective temperature between 
3550 and 3720K. 

From the tracks of D' Antona and Mazzitelli we get a mass 
of O.3OM and an age of 5 x 10 5 yrs for NX Pup C (see Table 
6 and Fig. 5). 

Our results are sensitive for the choice of evolutionary 
tracks. Tracks based on convection models with mixing-length 
theory yield a somewhat higher age for NX Pup C (10 6 yrs). 

4. Summary 

We investigated the properties of stars located close to the Her- 
big Ae/Be star NX Pup A. We found that its close companion, 
NX Pup B, is very likely a pre-main sequence star exhibiting 
strong IR excess, though the lack of a spectrum does not allow 
us to assess its mass and age. Nevertheless, its somewhat lower 
luminosity than NX Pup A, and its position in colour-colour- 
diagrams suggest a spectral type between mid F and late G. 



Table 6. Evolutionary status of NX Pup A, B, and C 



NX Pup 



B 



sep. 
PA 
SpT 
L/L Q 
age 



0'.' 128 ±0'.' 008 6'.' 98 ± 0'.'04 

62?4±5?7 45?3±0?2 

F5-G8 M0.5-M1.5 
9-18 0.45-0.85 

,6 



A7-F2" 
15-29 

5 x 10 6 yrs 0.3-5 x 10" yrs 5 x 10 b yrs 



a Brand et al. 1983, Reipurth 1983, Blondel & Tjin A 
Djie 1994 



The more distant companion, NX Pup C, which may or may 
not be physically associated to NX Pup A/B, is identified as 
a low-mass classical T Tauri star whose age around 5 x 10 5 
years is somewhat younger than the ages deduced by others 
for NX Pup AB and CGI. This provides the first evidence that 
low-mass star formation occurred in this cometary globule. 
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